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Abstract 

We compute the magnetic structure factor, the singlet correlation function 

and the momentum distribution of the one-dimensional Kondo lattice model 

at the density p = 0.7. The density matrix-renormalization group method 

is used. We show that in the weak-coupling regime, the ground state is 

paramagnetic. We argue that a Luttinger liquid description of the model 

in this region is consistent with our calculations . In the strong-coupling 

regime, the ground state becomes ferromagnetic. The conduction electrons 

show a spinless-fermion like behavior. 
PACS: 75.30.Mb, 75.30.Cr, 75.40.Mg 
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The heavy-fermion materials display a variety of low temperature-statesS. In some com- 
pounds, the Fermi-liquid state with huge quasiparticle mass is stable down to the lowest 
attainable temperatures. In others, the Fermi-liquid becomes magnetic or superconductor 
or both. The magnetic structures vary from simple Neel antiferromagnetism through in- 
commensurate order to ferromagnetism. This rich phenomenology arises from the interplay 
between a wide conduction-electron <i-band and a partially filled narrow /-band of rare-earth 
actinide or lanthanide elements. One of the canonical models for the description of these 
systems is the Kondo lattice model (KLM). 

After almost two decades of intense studies, the KLM is far from being completly 
understood!!. Even for the simplest case of one dimension, a consensus is reached only 
for the half-filled easel. The model is an insulator. A spin gap and a charge gap open for 
all non-zero values of the Kondo coupling. In the metallic phase, Monte carlo! and exact 
diagonalizationi calculations have found that the KLM presents two phases. A paramag- 
netic (PM) state with RKKY correlations is stable at low values of the Kondo coupling. 
The PM state evolves to a ferromagnetic (FM) ground state when the coupling becomes 
greater than a critical value. This FM phase is found to be numerically in accord with the 
results of strong-coupling expansion!. At very low densities, the model has a FM phase at 
all non-zero couplings^. The Monte Carlo study was, however, performed at finite tempera- 
ture while the exact diagonalization one was restricted to a lattice of 8 sites. Therefore, the 
trends displayed in these computations remain to be confirmed in the ground state of much 
longer chains . Using standard notations the KLM may be written as follows: 

H = -tJ2(4 s c i+ls + h.c.) + Ja'E S - • S i/ (!) 

is i 

where 

S? c = 2 E C ts a ss lC is'i &Tf = 2 E fis a< ss'fis' (2) 



s,s' 



with the constraint that 

E tifis = 1 (3) 



the letters c and / stand for conduction and localized electrons respectively. We applied the 
density-matrix renormalization group (DMRG) technique to this one-dimensional Kondo 
lattice. We chose a density p = 0.7, typical of those found in 3D compounds . We believe 
that the properties displayed at this density will reflect the behavior of the model in the 
whole region of moderate doping. Although the physical range of Jk corresponds to small 
couplings, it is also interesting to study the strong-coupling regime in order to understand 
the behavior of the model in the whole range of parameters. We varied Jk from 0.25 to 10. 
In the DMRG, an iteration of the algorithm consists in adding two sites at each step. It 
can be immediately realized that there is a problem in keeping the electron density fixed 
during the iteration process. To get around such a problem, we constructed the reduced 
density matrix from the two states whose electron numbers bracket the desired density. In 
the DMRG method, the states are also labelled by the ^-component of the total spin. In 
the present study, we work in the subspaces having = and ±|. For more details, the 
reader is referred to a recent paper by Chen and one of usi, where the method has been 
succesfully checked for the one-dimensional t — J model. The maximum lattice size we have 
reached is 75 with up to 180 states kept in the two external blocks. The truncation error 
of the Hilbert space at each iteration is around 7.10 -4 at Jk = 0.5 and less than 10 -6 at 
Jk = 10. We have computed the binding energy, the magnetic structure factor, the on- 
site conduction electron-localized spin correlation and the electron momentum distribution. 
Our results confirm the conclusions of Monte Carlo§ and exact diagonalizatiorJ! studies. A 
paramagnetic state is stable in the small coupling regime. This state is characterized by 
a maximum in magnetic structure factor at 2k p c = vrp and a singularity in the electron 
momentum distribution function at kp c . In the strong-coupling regime, the ground state is 
FM. The singularity of the conduction-electron momentum distribution is shifted to 2fcp c . 

When calculating the ground-state energy, we have taken the average of the lowest states 
having p\ and p2 such that p\ < p < p 2 (see Ref.9). This way, the ground-state energy per 
site of the non-interacting system can be reproduced up to four digits. In the Jk = oo 
case, the N c conduction electrons form perfect on-site spin-singlets with the localized spins. 
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The other N — N c /-spins remain free. The system presents a 2^ N ~ Nc ^ fold spin degeneracy. 
When the coupling is strong but finite, the ground-state energy per site is very close to 
— \JkP- The binding energy is defined as E B = [Eq{Jk = 0, N) — E G (J K , N)} /N, where 
Eg(Jk,N) is the ground-state energy. At a given density p, E B per site will be very close 
to -Eoo = fpJft- + e . The quantity e = — ^sin(|p) is the energy per site of the non- 
interacting case. In Fig. I, we show that our results are consistent with such an analysis. 
The convergence to E^ is very smooth. However Fig.l also indicates that this picture 
breaks down around Jjc = 2. Below this value of the coupling, the system enters in the 
small coupling regime characterized by pfJk < 1, Pf being the density of state at the Fermi 
level of the non interacting Hamiltonian. For the one-impurity Kondo problem, the binding 
energy is given by E\ = 3t(pp.JK) 2 Ln(2t) + At exp — j^ Ei. The non analytic part of E\ 
defines the Kondo temperature Tk- In the inset of Fig.l, we have compared E B with E\. 
One can see that E B is greater than E\. This enhancement of the binding energy of the 
lattice over the one-impurity case results from the intersite magnetic interaction. However, 
the latter conclusion does not necessarily mean that the Kondo temperature of the lattice 
is greater than Tk- 

The magnetic properties of the KLM are studied by calculating the spin structure factors 
of the conduction electrons, S c {k) and of localized spins, Sf{k). These quantities are defined 
here as follows: 

s cj(k) = jtJ2 ( s i, c ,/- s m,c,/) exp[i(7 - m)k\. (4) 

l,m 

At the first steps of the algorithm, pi and P2 are significantly different from p and boundary 
effects are non negligible. Thus, we started the calculation of the correlation functions when 
the lattice size was around 35. We have noticed that the value of the correlation functions 
are extremely sensitive to that of Sj> . We have not taken the average as for the ground- 
state energy because when the lowest state corresponding to pi has Sf. = 0, that of p2 is 
= ±|, and vice versa. We have used only the state with = 0. In order to reduce 
the effects due to the variation of the density during the iterations, we have started the 
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calculation of the correlation functions at different lattice size around 35 and then taken 
the average. In agreement with previous studies, we have found that the structure factor 
shows the competition of PM and FM phases. For small Jk, both Sf(k) (Fig. 2) and S c (k) 
(Fig.?>) have a maximum at 2kp c . The strong spin correlation observed in Sf(k) is due 
to the RKKY interaction. As the coupling is increased, this maximum flattens out. It 
completly disappears around Jk = 2. We note that this result compares well with that of 
exact diagonalization and Monte Carlo studies. When Jk is greater than this value, a new 
maximum arises at k = 0. It has been shown that in the strong-coupling regime, the effective 
interaction between the /-spins is ferromagnetic!!. The maximum at k = of S c (k) can be 
seen as a consequence of small ferromagnetic correlations tracking those of the /-spins. 

The suppression of the RKKY correlations can be understood by the Kondo mechanism. 
When Jk is small, the on-site singlet correlation (Sj^.Sjj) shown in Fig A, is different from 
that of perfect on-site singlets —\p- This is due to the non-local character of the Kondo 
singlets. This means that the singlet formed by an impurity and the conduction electrons has 
a spatial extension. Consequently, the RKKY interaction between impurities is favoured. 
As the coupling is increased, the size of the singlet reduces. It rapidly approaches that of 
perfect on-site singlets. As a result, the RKKY mechanism is suppressed. However contrary 
to the conventional view of the a global singlet state is unstable against a FM state 

in the strong Jk limit. In this limit, the local Kondo singlets are mobile. It is the undelying 
electron motion in reduced dimension that is responsible of the FM correlations. 

Now, we turn to the discussion of the conduction electron momentum distribution n(k). 
It has been argued that in the weak-coupling regime, the KLM may present a PM Lut- 
tinger liquid (LL) structure!. The LL behavior is characteristic of many one-dimensional 
interacting electron systems^. In the LL theory of PM systems, S(k) has a maximum at 
2kp c and n(k) presents a singularity at kp c . The electron momentum distribution has the 
form that follows: 

n(k) = rik F — C \k — kp\ a sgn(k — kp) (5) 
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We cannot calculate accurately the value of the exponent a. But Fig.h clearly shows that 
the singularity is located at k — kp c in the weak-coupling regime. The existence of this 
singularity is consistent with the presence of the maximum of S c j(k). So, the description of 
the model in the weak-coupling regime in terms of LL is plausible. Further investigations 
of the excitation spectrum are necessary to get the full answer. At intermediate Jk, n(k) is 
found to be very smooth so that it becomes hard to identify any singularity. At strong Jk 
however, the singularity as in a spinless fermion system. We believe 

that this is due to the action of the Kondo coupling which supress a double occupancy of the 
conduction electrons. It freezes the electron-spin degrees of freedom. The transition from a 
PM to a FM phase does not necessary mean that the LL description breaks down. It can 
be yet interpreted as a transition from a PM LL to a ferromagnetic LL. 

Finally, we touch upon the question of the size of the Fermi surface in the PM phase. In 
the weak-coupling region, the KLM is an effective model of the periodic Anderson model 
(PAM). In the PAM the / electrons are mixed to the conduction electron through an 
hybridization term. The PAM is believed to have a large Fermi surface containing both 
conduction electrons and localized electrons. In the KLM, however, there is not any hy- 
bridization between the two kind of particles. It is still a matter of debate wether or not the 
KLM has a large Fermi surface. A large Fermi surface supposes the existence of a maximum 
at 2k p c + 7r in the structure factor or of a singularity in the momentum distribution function 
at kp c + ~. Our results do not show any significant feature at this wave number. There is yet 
the possibility that the singularity at the position of the large Fermi surface is very smalllil. 

In conclusion, we have used the DMRG to study the KLM at the density p = 0.7. We 
believe that the behavior of the model at this density is characteristic of the moderate doping 
region. In agreement with previous exact diagonalization and Monte Carlo calculations, we 
have shown that the model presents a transition around Jk = 2. The weak-coupling region 
is PM. The magnetic structure factor has a maximum at 1kp c . The electron momentum 
distribution function displays a singularity at k Fc . We have argued that this is consistent 
with the description of the model in terms of LL. The strong-coupling region is FM. In 
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this phase, the singularity in the electron momentum distribution is shifted to 2kF c - This 
can be interpreted as the presence of a ferromagnetic LL structure. 
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grant from the Natural Sciences and Engineering research Council (NSERC) of Canada and 
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FIGURES 

FIG. 1. The binding energy of the Kondo lattice Eb (circles) versus the Kondo coupling. The 
dashed line corresponds to the strong-coupling limit. In the inset Eb is compared with the binding 
energy of the one-impurity (diamonds) Kondo problem. 

FIG. 2. The magnetic structure factor of the localized spins for p = 0.7 at Jk = 0.5 (circles), 
Jk = 1 (diamonds), Jk = 2 (stars) and Jk = 4 (triangles). 

FIG. 3. The magnetic structure factor of the conduction electrons for p = 0.7 at Jk = 0.5 
(circles), Jk = 1 (diamonds), Jk = 2 (stars) and Jk = 4 (triangles). 

FIG. 4. The on-site conduction electron-localized spin correlation versus the Kondo coupling. 
The dashed line corresponds to the strong-coupling limit 

FIG. 5. The electron momentum distribution n(k) for p = 0.7 at Jk = 0.5 (circles), Jk = 1 
(diamonds), Jk = 2 (stars) and Jk = 4 (triangles). 
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